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ABSTRACT 



We use numerical simulations of isolated galaxies to study the effects of stellar feedback on the for- 
mation and evolution of giant star-forming gas "clumps" in high-redshift, gas-rich galaxies. Such galactic 
disks are unstable to the formation of bound gas-rich clumps whose properties initially depend only on 
global disk properties, not the microphysics of feedback. In simulations without stellar feedback, clumps 
turn an order-unity fraction of their mass into stars and sink to the center, forming a large bulge and kick- 
ing most of the stars out into a much more extended stellar envelope. By contrast, strong radiative stellar 
feedback disrupts even the most massive clumps after they turn ~ 10 — 20% of their mass into stars, in a 
timescale of ~ 10 — 100 Myr, ejecting some material into a super- wind and recycling the rest of the gas 
into the diffuse ISM. This suppresses the bulge formation rate by direct "clump coalescence" by a factor 
of several. However, the galactic disks do undergo significant internal evolution in the absence of merg- 
ers: clumps form and disrupt continuously and torque gas to the galactic center. The resulting evolution is 
qualitatively similar to bar/spiral evolution in simulations with a more homogeneous ISM. 
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1 INTRODUCTION 



In recent years, considerable attention has been paid to what ap- 
pears to be a population of massive (M* ^ 10^^ — 10^^ Mq), gas- 
rich (Mgas/(Mgas +M*) ^ 0.3 — 0.7), rapidly star forming (M* ^ 
10 — 200 Mq yr~^) "disk" galaxies (flattened and apparently rota- 
tionally supported, albeit with the ratio of rotational to random ve- 
locity V/a r^a. few) at redshifts z 2 — 3 (see e.g. "Chapman et alj 
[20051 iGreveeTaT 2005, Tacconi et al. 2006, Genzel et al. 20^ 
Tacconi et al. 2010, Shapiro et al. 2009, Da ddi et al.|20T0|[F5rster| 
jSchreiber et al.,2006^ 2011b) . The most extreme star-forming pop- 
ulations at these redshifts (with M* ~ 100 — 3000 Mq yr~^) appear 
to be major merger-induced starbursts, but these "disky" systems 
dominate the intermediate-luminosity population which includes 
a large fraction of the total SFR density (see references above 
and [Tacconi et al.|2008j ). The morphology of the young stars (ob- 
served in the rest-frame UV and optical light) and star-forming or 
molecular gas (e.g. Ha) is characteristically irregular ("clumpy" or 
"clump-chain" morphologies), with a significant fraction (tens of 
percent) of the light in ~a few massive (^ 10^ Mq), ~kpc-scale 
clumps or blobs ( Griffiths et al.'1994',"Cowie et al."1995»,'Giavahsco' 
etal . 1996, Abraham et al. 1996, Elmegreen et al. 2004b a 2005b, 
Kriek et al. 2009 , Genzel et aL|2011| |Swinbank et al.|2011| . 



Theoretically, the origin of these morphologies is controver- 
sial. Some early models argued they were irregular because of on- 
going mergers ( [Somerville et al.|2001| ); observations of local ana- 
logues (and studies in which these systems are mock-observed as if 
they were at high redshifts) suggest that this is true for least some 
fraction of the observed systems ( [Gverzier et al.|2010]|2009| [Petty] 



[et al.|[2009] ). However, the constraints on the structure and kine- 
matics of the massive z ^ 2 systems suggest that a sizeable frac- 



matics 01 tne massive z ^ z systems sugg i 
tion are more quiescen{^ |Shapiro et al. 2 



2008 



Forster Schreiber 



[et al.|2006[|2011aj, although this is n ot definitive evidence against 
mergers ( [Robertson & Bullock|2008[[Hammer et al.[2009| ). And the 
merger rate may not be sufficiently high to explain the abundance 
of these intermediate-SFR systems (e.g. Dekel et al. 2009b , Stewart] 
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[etaLl2 009 Hop kins et al.|2010[[HopEns & Hernquist 2010}. 

An alternative explanation for the observed morphologies is 
that dumpiness stems from disk fragmentation. Any gas disk in 
which cooling is efficient (cooling time less than the dynami- 
cal time so that turbulent support dissipates in a crossing time, 
which is easily satisfied here) will fragment at the Toomre scale, 
Rt - a^/ivGT. or Mr - a^TrG^T. ( ,Toomre||1964| ); for disks 
which equilibrate at marginal stability (Toomre g ~ 1), this is just 
Rt ^ /gas^disk (Mr ^ /gasMdisk)[^In the Milky Way, these clumps 
correspond to massive GMCs (R ^ lOOpc, M ^ IO^Mq). But in 
a massive, gas-rich disk, the scale-length and mass of the resulting 
clumps can be as large as ^kpc and ^ 10^ — IO^Mq, respectively, 
similar to the structures observed. 

This instability has been well-known for ^ 40 years. The 
implications of such clumping in the ISM are less clear, how- 
ever. If there were no feedback from stars to expel mass from 
GMCs/clumps, they would (locally) collapse in a few free-fall 



^ By "quiescent" here, we mean non-interacting/non-merging, as opposed 
to "quenched" or non-star forming. 

^ We have dropped numerical pre-factors that depend on the detailed disk 
structure. Also we could take /coid instead of /gas to reflect the dynamically 
thin/cold stellar components, although in detail the scaling for a gas plus 
stellar system, while qualitatively similar, becomes much more complex. 
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times and turn most of their mass into stars. They would then sink 
to the galaxy center under the influence of dynamical friction on a 
timescale ^sink ^orbit Mdisk/Mr; for GMCs, this would take a Hub- 
ble time, but for massive clumps in the high redshift systems of 
interest, this could occur in a few disk orbital periods 30 — 60 
clump dynamical times, or ^ 0.5 Gyr). This sinking would quickly 
build up a ^kpc-scale bulge ( Shlosman & Noguchi|1993[|Noguchi| 
\\999) . This is also well-known and has been seen since the first 
simulations of very gas-rich disks ( Hernquist ]1989 [Barnes & Hern-| 
|quist'1991;"Shlosman & Noguchi 1993 , Barnes & Hernquist|1996| 
Sommer-Larsen et al. 1999, Robertson et al. 2004]l 

However, in the MW and nearby galaxies, typical GMCs are 
disrupted by stellar feedback after a few clump dynamical times 
(a few Myr), and recycle ~ 95% of their mass back into the ISM 
([Zuckerman & Ev ans] 1974] [Williams & McKee| 19971 [Evans 1 1999) 
[Evans et a l. 2009). This appears to be true for both low and high 
mass clouds, including local examples with masses and sizes very 
similar to those observed at high redshifts (Mdoud ^ 10^ — IO^Mq; 
[Wilson et al.|2003[[2006[ and references therein). 

If Toomre mass clumps were dispersed after a few dynami- 
cal times, it may be more accurate, in simulations, to treat the ISM 
as a "smooth" medium (i.e. not resolve the formation of individ- 
ual clumps) rather than resolve clump formation but not their de- 
struction/recycling. In fact, the one-way fragmentation seen in sim- 
ulations without feedback was one of the primary motivations for 
modelers to include stronger prescriptions for stellar feedback in 
numerical galaxy models. Lacking the ability to explicitly resolve 
the processes that destroy clouds, this commonly amounted to as- 
signing the gas a much larger "effective pressure" than its nominal 
thermal pressure so that runaway collapse could not occur. Such an 
effective pressure was intended to mimic a global average turbu- 
lent pressure maintained by a sub-grid process of continuous clump 
formation and destruction f Hemquist|1989[ [Springel & Hernquist[ 
[2003 ). This was also motivated by cosmological simulations - the 
one-way break-up of gas-rich disks was a significant factor pre- 
venting the formation of realistic disk galaxies and led to simula- 
tions which uniformly predicted bulge-to-disk ratios significantly 
larger than observed ("Sommer-La rsen et al.|1999[ [Rob ertson et al." 
[2004[ [Governato et al. 2004 2007)). With the recent observations 
of these "clumpy" high-redshift galaxies, however, renewed inter- 
est has been paid to "weak feedback" models of the I SM pmmeli 



etal. 2004, Bour naudet al.[2007|[B ournau d erar.[2008[[Elmegreen 
et al.|2008[[Agertz et al.|2009|[Dekel et al.|2009a|[Ceverino et al. 
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Without an a priori resolved model of stellar feedback, the 
dominant feedback physics and their effects remain uncertain. 
[Genel et ar] |2012) considered the evolution of gas-rich disks in- 
cluding various phenomenological models in which stellar feed- 
back was parameterized by a "wind efficiency" (in short, for ev- 
ery mass AM* in stars formed, a mass 77 AM* is automatically 
kicked out of the galaxy). They showed that for order-unity wind 
efficiencies of the sort needed in cosmological models to explain 
the shape of the galaxy mass function ( Op penheimer et al.|2010| ), 
clumps form but are disrupted by feedback before turning a large 
fraction of their mass into stars and so do not spiral to the center. 
Observations have now suggested comparable wind mass-loading 
factors in the clumps ( Genzel et al.[[201l| [Newman et al.[[20T2| . 
These models, while a powerful "proof of concept," are not able 
to predict the wind properties or detailed ISM structure. The stud- 
ies of clumpy disks that have included physically-motivated models 
of stellar feedback have focused only on thermal heating from SNe; 
however, at the high gas densities of typical "clumps" (> 100 cm~^) 



gas which is shock-heated by SNe ejecta will cool in a timescale 
> 100 times shorter than its dynamical time, so these models un- 
surprisingly find similar results (in terms of clump survival) to 
feedback-free cases ( Elmegreen et al.|2008[[Ceverino et al.|2010| ). 
And at higher densities, in fact, the SNe timescale becomes longer 
than the local collapse/free fall time, so they cannot act efficiently 
(Evans et al.|200 9 ). 

A large number of feedback mechanisms may drive turbulence 
in the ISM and help disrupt GMCs, including: photo-ionization, 
stellar winds, radiation pressure from UV and IR photons, proto- 
stellar jets, cosmic rays, supernovae, and gravitational cascades 
from large scales (e.g. Mac Low & Klessen 2004 and references 
therein). In ^Hopkins et al. (2011) (Paper I) and Hopkins et al.[ 
( [2012b| ) (Paper II) we developed a new set of numerical models 
to incorporate feedback on small scales in GMCs and star-forming 
regions, in simulations with pc-scale resolutionj^These simulations 
include the momentum imparted locally (on sub-GMC scales) from 
stellar radiation pressure, radiation pressure on larger scales via the 
light that escapes star-forming regions, HII photoionization heat- 
ing, as well as the heating, momentum deposition, and mass loss 
by SNe (Type-I and Type-II) and stellar winds (O star and AGB). 
The feedback is tied to the young stars, with the energetics and 
time-dependence taken directly from stellar evolution models. Our 
models also include realistic cooling to temperatures < lOOK, and 
a treatment of the molecular/atomic transition in gas and its effect 
on star formation. We showed in Papers I & II that these feedback 
mechanisms produce a quasi- steady ISM in which giant molecular 
clouds form and disperse rapidly, after turning just a few percent 
of their mass into stars. This leads to an ISM with phase structure, 
turbulent velocity dispersions, scale heights, and GMC properties 
(mass functions, sizes, scaling laws) in reasonable agreement with 
observations. In [Hopkins et al.[p012a| ) (Paper III), we showed that 
these same models of stellar feedback produce the elusive winds 
invoked in almost all galaxy formation models; the combination 
of multiple feedback mechanisms give rise to massive, multi-phase 
winds having a broad distribution of velocities, with material both 
stirred in local fountains and unbound from the disk. 

In this paper, we use these models to study massive, GMC- 
analogue clumps in high-redshift star-forming disks, and their im- 
plications for bulge formation. 

2 THE SIMULATIONS 

The simulations used here are described in detail in Paper I (Sec. 2 
& Tables 1-3) and Paper II (Sec. 2). We briefly summarize the most 
important properties here. The simulations were performed with the 
parallel TreeSPH code GADGET-3 ( [Springel[[2n05] ). They include 
stars, dark matter, and gas, with cooling, star formation, and stellar 
feedback. 

Gas follows an atomic cooling curve with additional fine- 
structure cooling to < 100 K, with no "cooling floor" imposed. Star 
formation is allowed only in dense regions above n > lOOOcm"^, 
at a rate p* = e pmoi/fe where is the free-fall time, pmoi = /h2 P is 
the molecular gas density, and e = 1.5% is the observed efficiency 
at these densities ( [Krumholz & Tan|2007l ). We follow [Krumholz &[ 
[Gnedin[(2011| ) to calculate the molecular fraction in dense gas 
as a function of local column density and metallicity. In Paper I and 
Paper II we show that the galaxy structure and SFR are basically 



^ Movies of these simulations are available at https : / /www . c f a . [ 
harvard . edu/~phopkins/ Site /Re search . html 
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Figure 1. Morphology of the gas and stars in a simulation of a massive z ~ 2 — 4 rapidly star-forming disk with M* ~ 100 Mq yr~^ . Images are shown at ~ 2 
orbital times when the disk is in a feedback-regulated steady-state. We show face-on (upper) and edge-on (lower) projections. We compare two simulations: 
our "standard" model with all feedback mechanisms enabled (left two columns), and a model with identical conditions but radiation pressure (the dominant 
feedback mechanism) disabled {right two columns). Top: Stars. The image is a mock ugr (SDSS-band) composite, with the spectrum of all stars calculated 
from their known age and metallicity, and dust extinction/reddening accounted for from the line-of-sight dust mass0For each model, the left column is an 
image with effectively infinite resolution and depth (brightness follows a logarithmic scale with a stretch of ~ 3dex); the right column is the same image, 
convolved with a Gaussian PSF typical of HST WFC3 observations (FWHM 1.4kpc = 0.16^^ at z = 2) and with a linear stretch. Bottom: Gas. Brightness 
encodes projected gas density; color encodes gas temperature with blue/white being T < 1000 K molecular gas, pink ~ 10^ — 10^ K warm ionized gas, and 
yellow > 10^ K hot gas. Again, left is at extremely high resolution/depth (logarithmically scaled with a ~ 6dex stretch from ~ 10~^ — 10^ Mq pc~^); right 
is the same with PSF convolution typical of e.g. SINFONI AO (FWHM 1.7kpc = 0.2^^ at z = 2) and linear stretch. Gravitational collapse forms massive 
kpc-scale star cluster complexes that give rise to the clumpy morphology (edge on, similar to "clump-chain" galaxies) that is further enhanced in the optical by 
patchy extinction. The clumps have considerable sub-structure. With all feedback included, strong outflows are present, emerging both from the complexes and 
the disk as a whole, driven by the massive star formation rate. These outflows maintain a thick gas disk and disrupt many of the clumps, which continuously 
re-form. With radiation pressure disabled, the outflows are much weaker (some "venting" hot gas), the clumps collapse further, are more clearly rotationally 
supported, and leave less of a diffuse low-brightness stellar background (having coalesced more efficiently). 



independent of the small-scale SF law, density threshold (provided 
it is high), and treatment of molecular chemistry. 

We consider an initial galaxy model intended to represent an 
(intentionally extreme) high-redshift massive "starburst disk" (typ- 
ical of high-luminosity BzK galaxies or low-luminosity SMGs - 
intermediate SFR systems - at z ~ 2 — 4). The disk has a small 
{B/T < 0.05) initial H emquist] ( |1990| ) profile bulge (M^ = 5 x 
10^ Mq (chosen to be small to avoid suppressing disk instabilities 
and clearly isolate new bulge growth), scale-length a = 1.2kpc), an 
exponential stellar disk (Md = 3 x 10^^ Mq \rd — \ .6kpc with sech^ 
profile scale-height h — 130pc) a gas disk (Mg = 7 X 10^^ Mq; 
rg = 3.2kpc), and a dark matter halo (Mhaio = 1.4 x 10^^ Mq, con- 



centration c = 3.5, scaled to lie on the mass-concentration relation 
for z = 2 halos). The model has ^ total particles and 5pc grav- 
itational softening. 

Stellar feedback is included, via a variety of mechanisms. 

(1) Local Momentum Flux from Radiation Pressure, Super- 
novae, & Stellar Winds: Gas within a GMC (identified with an on- 
the-fly friends-of-friends algorithm) receives a direct momentum 
flux from the stars in that cluster/clump. The momentum flux is 
P = ^SNe + + ^rad, whcrc the separate terms represent the direct 
momentum flux of SNe ejecta, stellar winds, and radiation pres- 
sure. The first two are directly tabulated for a single stellar popula- 
tion as a function of age and metallicity Z and the flux is directed 



© 0000 RAS, MNRAS 000, 000-000 



4 Hopkins et al. 




Figure 2. Statistics of gaseous clump lifetimes and star formation efficiencies, in simulations with and without stellar feedback. We identify all gas clumps that 
can be well-resolved (> 100 particles) down to 10% of the clump mass (> 1000 particles at peak); lifetime is defined as the time above 10% the maximum 
cloud gas mass. Left: Integrated star formation efficiency (stellar mass formed in clump over maximum cloud gas mass). We compare the "standard" model 
(all feedback included) to a model with no feedback. Middle: Cloud lifetime in units of the (mass-weighted average) free-fall time fff ~ 0.54/ v^Gp within 
the clump. Right: Lifetime versus maximum cloud gas mass. We show the result for individual clouds (small black and large red points for the standard/no 
feedback models, respectively) and the lines show the ±1 cr range at each cloud mass for each model. With feedback present, even the most massive gas clouds 
live for a few dynamical times and turn ~ 10% of their mass into stars. Without feedback, clouds persist until most of their dense gas is turned into stars, 
giving star formation efficiencies of ~ 50%; their lifetimes are longer and they "end" only because the gas is exhausted (the relic stellar clumps are much 
longer-lived). 
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Figure 3. Remnant stellar mass profiles for our simulations, with differ- 
ent stellar feedback mechanisms enabled or disabled (as labeled). The pro- 
files are averaged over ~ 100 projections. The initial disk profile is shown 
for comparison. There is internal (rapid "secular") evolution in all cases: 
gas that loses angular momentum builds a pseudo-bulge while some gas 
gains angular momentum, making the outer disk larger. With our standard 
feedback model, the bulge is small (B/T ~ 0.1) and disky (jis ~ 1), and 
the outer disk does not evolve significantly. Removing SNe heating makes 
little difference. Removing radiation pressure, however, is nearly identical 
to removing all feedback, which yields a much larger, more concentrated 
{us ~ 2) bulge as well as a stellar envelope at large radii that includes 
~ 60 -70% of the mass. 



away from the stellar center. Because this is interior to clouds, the 
systems are always optically thick, so the latter is approximately 

A^d ~ (1 + TiR)Lincident/c, whcrC 1 + TiR = 1 + Sgas aCCOUntS 

for the absorption of the initial UV/optical flux and multiple scat- 
terings of the IR flux if the region is optically thick in the IR (with 
Egas calculated for each particle given its location in the clump). 



(2) Supernova Shock-Heating: Gas shocked by supernovae 
can be heated to high temperatures. We tabulate the SNe Type-I and 
Type-II rates from [Mannucci et~aLl ( [20061 ) and STARBURST99, 
respectively, as a function of age and metallicity for all star particles 
and stochastically determine at each timestep if a SNe occurs. If so, 
the appropriate mechanical luminosity is injected as thermal energy 
in the gas within a smoothing length (nearest 32 gas neighbors) of 
the star particle. 

(3) Gas Recycling and Shock-Heating in Stellar Winds: 

Gas mass is returned to the ISM from stellar evolution, at a rate tab- 
ulated from SNe and stellar mass loss (integrated fraction 0.3). 
The SNe heating is described above. Similarly, stellar winds are 
assumed to shock locally and inject the appropriate tabulated me- 
chanical luminosity L(^, Z) as a function of age and metallicity into 
the gas within a smoothing length. 

(4) Photo-Heating of HII Regions: We also tabulate the rate 
of production of ionizing photons for each star particle; moving ra- 
dially outwards from the star, we then ionize each neutral gas par- 
ticle (using its density and state to determine the necessary photon 
number) until the photon budget is exhausted. Ionized gas is main- 
tained at a minimum ^ 10^ K until it falls outside an HII region. 

(5) Long-Range Radiation Pressure: Photons which escape 
the local GMC (not absorbed in mechanism (1) above) can be 
absorbed at larger radii. Knowing the intrinsic SED of each star 
particle, we attenuate integrating the local gas density and gradi- 
ents to convergence. The resulting "escaped" SED gives a flux that 
propagates to large distances, and can be treated in the same man- 
ner as the gravity tree to give the local net incident flux on a gas 
particle. The local absorption is then calculated integrating over a 
frequency-dependent opacity that scales with metallicity, and the 
radiation pressure force is imparted (and luminosity removed). 

In implementing (l)-(5), all energy, mass, and momentum- 
injection rates are taken from stellar population models ( Leitherer| 
|et al.|1999] ), assuming a ,Kroupa ( 2002) IMF, without any free pa- 
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rameters. More details, numerical tests, and resolution studies (up 
to 10^ particles with 3.5 pc softening lengths) for these models are 
discussed in Paper 11. 

We note that some recent studies of low-resolution cosmo- 
logical simulations comparing GADGET and the moving mesh 
code AREPO { |Springel||2010| ) have highlighted some differences 
between smoothed particle hydrodynamics and grid methods for 
some cosmological inflow problems ( |Vogelsberger et ar]|2011[ 
[Keres et al.|2012[[Torrey et al.|2011| ). However, we have also per- 
formed idealized simulation comparisons between the individual, 
high-resolution galaxy models here (as well as galaxy mergers) 
and found excellent agreement for e.g. the star-formation rates, gas 
inflow through the disk, and surface density profiles; these com- 
parisons will be presented in |Hayward et al.|f2011| ). We have also 
re-run the simulations here with an alternative density-independent 
formulation of SPH (which resolves these discrepancies) as pro- 
posed in Sait oh & Maki no (2012); again we find excellent agree- 
ment. This owes in part to the resolution, as well as to the fact that 
the code differences are minimized when flows are super- sonic, al- 
ways true here (see |Bauer & Springel|2012| ). To the extent that there 
are differences, the moving-mesh methods may tend to slightly 
shorter clump lifetimes because stripping and mixing of the gas 
as it moves are more efficient (and artificial shock dissipation is 
reduced) (see |Sijacki et al.|2012| ). 

3 RESULTS 

Fig. [l] compares the gas and stellaij^ morphologies of runs with 
varying feedback properties. In all cases, the disks form massive 
(> IO^Mq) ^kpc-scale clumps, which in turn form massive star 
cluster complexes. The gas densities are sufficiently high that ex- 
tinction can lead to a clumpy optical morphology as well (by 
entirely extincting out regions in the galaxy). The dumpiness is 
prominent in the gas and stars (although, with feedback present, 
there is a low surface brightness smooth stellar disk that is not 
visible at the present observational depth). Seen edge-on, the gas 
can resemble a "clump-chain" morphology. With feedback present, 
strong outflows arise from throughout the disk, driven by a massive 
star formation rate of > 1 00 Mq yr~^ 

In previous papers, we present detailed comparisons between 
our simulated disk properties and observations. The model disk 
masses, rotation curves, scale lengths, and gas fractions (see § [2| 
are chosen to match those observed. Their resulting SFRs (50 — 
3OOM0yr~^) are consistent with the Kennicutt- Schmidt relation 
(Paper I; Fig. 11 and Paper II; Fig. 7). Moreover, the simulated 
galaxies self-regulate at Toomre 2 ~ 1, with typical line-of-sight 
gas velocity dispersions a 50 — 80kms~^ and star-forming gas 
disk scale heights \ kpc or h/R ^ 0.2 (Paper II; Fig. 9). How- 
ever, we show in Paper II that our model disks always reach Q'-^ \ 
and so - for fixed global disk properties - a is almost completely 
independent of the SFR density and the microphysics of feedback. 
The maximum clump masses are ^ 10^ Mq with a few clumps per 
disk within a factor of a few of this mass (Paper II; Fig. 13); the 



corresponding clump sizes are ^ 0.3 — 2 kpc (Paper II; Fig. 15) giv- 
ing typical clump surface densities 300 — 1 000 Mq pc~^ (Paper 
II; Fig. 14). The clumps are marginally bound (and not rotation- 
ally supported, although mildly anisotropic), giving them internal 
velocity dispersions nearly equal to that of the "background" disk 
gas (Paper II; Fig. 16; see also | Hopkins|2012| ). The winds driven 
off of the disk and clumps have a broad velocity distribution from 
^ 200 — 750 kms~^ (Paper III; Fig. 4), with globally averaged out- 
flow rates in the high-speed superwind of Mwind ~ 0.5 — 2M* (Pa- 
per III; Figs. 6-7). All of these properties agree reasonably well with 
those observed (see |Genzel et al.|2011| and references therein). As 
shown in Paper II and Paper III, the simulations quickly establish 
these properties in the first couple dynamical times and maintain 
them in quasi- steady state as long as they are evolved, until gas is 
exhausted. 

We now quantify the properties of the clumps in more de- 
tail and discuss their implications for the internal ("secular," but 
rapid) evolution of the galaxy. We identify gas clumps following 
Paper II; specifically, we apply the sub-halo finder SUBFIND to 
the gas, which employs a friends-of-friends linking algorithm with 
an iterative procedure to robustly identify overdensities (for details 
and tests, see |Springel et al.||2001| ). In Paper II we show that vi- 
sual inspection and calculation of the binding energy of clumps 
confirm that this correctly identifies the obvious clumps; changing 
linking lengths and other numerical quantities has a weak effect on 
the properties here. Because the gas particles are unique, we can 
link clumps in time between different narrowly spaced snapshots 
(A^ 10^ yr) to construct a clump "merger tree We then obtain 
the growth/decay in mass with time for each clump; we define the 
clump "lifetime" as the time spent above 10% of the maximum 
clump mass (this cut is arbitrary, but not critical to our results). We 
can also define the total clump star formation efficiency by integrat- 
ing the SFR of all clump- associated particles at each time. 

Fig. |2] plots the distribution of gas clump lifetimes and star 
formation efficiencies. As shown in Paper II, the most massive 
gas clumps do tend to be longer-lived; however, with feedback 
present their lifetimes are still just a few internal dynamical times 
(^ 10 — lOOMyr; this is consistent with the - highly uncertain 
- observed stellar populations in clouds in e.g. Elmegreen e t aL| 
2009 , Forst erSchreiber et al.|201 lb] but more detailed observations 
would present a useful constraint on the predictions here). This is 
short compared to their inspiral time by dynamical friction. About 
^ 5 — 20% of the peak gas mass of the clumps is turned into stars; 
this is a factor of a few higher than the star formation efficiency 
of local GMCs (as expected by theoretical arguments; see [Murray] 
|et al.pOlO] ), but much less than the order-unity fraction needed to 
conserve the "relic" stellar clump as a bound entity]^ 

Without feedback, the gas collapses to higher densities 
rapidly, so the free-fall time oc p~^^^ becomes very short and 
the gas lifetime in units of that free-fall time becomes as long as 
^ 10 — lOOfe. This runaway collapse leads to large star formation 
efficiencies - 20 - 100%, with a mean of - 50% (it is not 100% 
because there is always some low-density material associated with 



^ The stellar luminosity in each band is calculated from each star particle 
according to the STARBURST99 model given its age, mass, and metallic- 
ity (and smoothed over the appropriate kernel). We then attenuate the stars 
following the method of Hopkins et al. ( 2005): we calculate the total dust 
column (from the simulated gas) along the line-of-sight to each star parti- 
cle for the chosen viewing angle (assuming a constant dust-to-metals ratio, 
i.e. dust-to-gas equal to the MW value times Z/Zq), and apply a MW-like 
extinction and reddening curve (as tabulated in|Pei 1992J. 



^ Here the "main branch" is defined by the most massive clump in the tree 
at each time. Mergers of other clumps are generally small - most growth is 
via accretion of "non-clump" gas, so we neglect the "merged branches." 
^ Many of the lowest mass-clouds, in both the models with and without 
feedback, have short lifetimes. The reason is, as shown in Paper II, many 
of these are never bound but are simply transient overdensities. We also 
caution that they may still be subject to resolution effects. But in any case 
these represent a small fraction of the total gas mass and SFR. 
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clouds). Note that the Ufetime of the "gas clump," defined as we 
consider it, may then still be relatively short, ^ 20 — lOOMyr; this 
is because most of the mass has turned into stars, a large fraction of 
which remains self-gravitating and will survive in a bound stellar 
clump until reaching the galaxy center. 

Fig. |3] shows the resulting mass profiles at ^ = 500 Myr (for 
the initial disk model, %n = Re/Vc{Re) ~ 12Myr, so this is 40 (7) 
dynamical (orbital) times at Re). By this point the gas is largely 
exhausted or blown out in all models and the evolution has slowed 
down; most "sinking" clumps have coalesced. That is not to say this 
is a "final state" - the system undergoes continuous long-term colli- 
sionless stellar/halo internal and secular evolution, and in a cosmo- 
logical context it would evolve significantly over ^Gyr timescales. 
We are simply comparing the relative effects of including or not 
including feedback. 

In all the simulations, there is some internal/"secular" evolu- 
tion away from the initial disk conditions. The gas is still asym- 
metric and turbulent, hence there are gravitational torques (see e.g. 
[Gammie 200 l||Hopkins & Qua taert 20lT] |Dekel et al.|2009a||Bour^ 
|naud et al.|2011[ [Forbes et al.|201 2 ). Some mass is channeled to 
the central ~kpc, building up the "bulge"; the angular momentum 
loss required for this is absorbed by the outer disk (basically the 
equivalent of an outer Lindblad resonance), which then migrates 
outwards into a more extended, flattened disk "envelope" (which 
may, in some cases, be directly observed; [Elmegreen & Elmegreen| 
12006] ). 

With feedback, this evolution is modest; fitting a bulge-disk 
decomposition to the profile (with a Sersic-profile bulge and ex- 
ponential disk), we obtain (B/T, Us, Rd) = (0.10, 1.1, 4.3 kpc) for 
the bulge mass fraction, Sersic index (lower Us < 2 representing 
more disk-like "pseudobulges"), and disk effective radius (compare 
to the initial Rd = 2.1 kpc) |j With no feedback, we see a dramatic 
central rise in the profile indicative of a large bulge, and a much 
more pronounced outwards evolution of the disk; (B/T, Us, Rd) = 
(0.37, 1.9, lOkpc). Fig. [3] compares several intermediate cases as 
well: if we remove radiation pressure, including only sources of gas 
"heating" (SNe, stellar wind shock-heating, and HII heating), the 
evolution is only slightly less severe than the feedback-free case. 
But if we include only radiation pressure, the profile is only some- 
what more evolved than the standard feedback case. If, however, 
we remove long-range radiation pressure forces, keeping only lo- 
cal radiation pressure in regions optically thick in the IR, then it is 
again more similar to the feedback-free case. As we show in Paper I 
and Paper III, the long-range photon momentum is critical for driv- 
ing winds; the short-range pressure tends to drive local turbulence 
within clumps, but may not completely disperse them. 

4 DISCUSSION 

We have simulated the isolated internal (rapid secular) evolution 
of gas-rich (initial /gas ~ 0.6), massive (10^^ Mq), disk-dominated 

^ We caution that the best-fit bulge profile (both with and without feedback) 
is quite sensitive to the exact fitting procedure, since these galaxies (being 
clumpy and asymmetric) do not have a well-defined "center" at the sub-kpc 
level (so it makes a large difference whether we center on a clump density 
maximum or between clumps). Here we use an iterative elliptical aperture 
density smoothing to isolate the central density maximum (see jHopkinsj 
|et a l. 2009 1, which gives the maximum "bulge" component of any fitting 
methods we explore. Our qualitative conclusions on relative bulge forma- 
tion efficiencies are the same with each method. The outwards evolution of 
the disk, however, is a much more robust indication of the degree of internal 
evolution and is completely insensitive to these fitting details. 



galaxies, designed to resemble the most massive star-forming disks 
observed at redshifts z^2 — 3. Our simulations include explicit, 
physical models for stellar feedback from a variety of sources: ra- 
diation pressure (from UV and IR photons), HII photoionization, 
SNe, and stellar winds. In Paper II, we show that these models 
agree reasonably well with observed systems in a wide range of 
measurable properties: effective radius, mass, gas fraction, and ro- 
tation curve (by construction), as well as gas velocity dispersion, 
SFR, Toomre Q parameter, gas disk scale height (h/R^ 0.2), and 
characteristic clump masses and sizes. Here we have quantified the 
galaxy morphology and the evolution of its stellar mass profile with 
different feedback mechanisms included or removed from the sim- 
ulations. 

We find that the models all experience some degree of internal 
evolution. This is expected: they are all gas-rich, disk-dominated 
systems in which the cooling time is short relative to the dynami- 
cal time, so are both locally and globally gravitationally unstable. 
We showed in Paper II that, independent of the details of feedback, 
the disks converge to Q ^ 1, which means that the gas velocity 
dispersion in all cases is essentially identicalj^This means that the 
characteristic Toomre length and mass scales are set by global prop- 
erties and are similar in all of our simulations, even with very dif- 
ferent feedback physics: Mxoomre ^ h^^gas ~ /gas^disk ^ IO^Mq. 
The size and mass scales for the most massive clumps are therefore 
similar in all our simulations. In Paper II we explicitly show that the 
mass function of gas clumps is also broadly similar. These clumps, 
then, are the generic extension of GMCs in the MW (which has a 
lower Toomre mass) to more gas-rich massive systems (although 
with some subtle distinctions; see |Krumholz et al.|2012| ). 

However, the effects of the resulting clumps on bulge forma- 
tion depends strongly on the inclusion of feedback. With no feed- 
back, the clumps collapse rapidly and turn ^ 50% of their mass 
into stars. They have no alternative but to sink to the center in a 
few orbital periods: this gives rise to a massive, ^kpc-scale bulge, 
and causes substantial migration of the remaining disk gas and stars 
to larger radii (since the angular momentum of the sinking clumps 
must be taken up by the remaining gas/stars, and tends to do so 
near the equivalent of an outer Lindblad resonance). These effects 
are, broadly speaking, well-known results of strong internal/secular 
evolution, though some of the details differ in this very clumpy 
limit (and may differ further with ongoing accretion). For an iso- 
lated disk, this evolution is self-terminating once sufficient bulge 
mass is formed, since the outward disk migration and inner bulge 
formation both stabilize the system; this generically leaves a large 
low-density disk "envelope" rather than a pure bulge (a well-known 
result shown in |Elmegreen et al.pOOSallBournaud et al.|2007| ). 

In our standard feedback model - with all mechanisms in- 
cluded - even the most massive clumps tend to be disrupted after 
they turn ^ 5 — 20% of their mass into stars. Even if those stars 
can remain bound after the clump is dispersed (in which case they 
will still tend to sink to the center, albeit on a factor ^10 longer 
timescale) this suppresses the efficiency of bulge formation via 
clump coalescence by an order of magnitude. And indeed we see 
much less dramatic signatures of internal/secular evolution: both a 
smaller bulge, and much weaker evolution of the outer disk. The re- 



^ This does not mean the dispersion cannot be "driven by" feedback; it 
simply means that once systems reach Q > I, excess energy is lost (e.g. 
in winds) or dissipated, without new star formation, and once Q < I, they 
collapse until sufficient stars form (or until sufficient gas is exhausted) to 
bring Q back up to unity. 
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suiting evolution is less rapid compared to one-way fragmentation, 
and produces a more traditional pseudobulge. However, the clump 
morphologies are not radically different; if anything, the clumps 
appear somewhat more prominent with feedback, because they are 
prevented from internally collapsing to small sizes. They are also 
less rotationally supported in better agreement with the observa- 
tions fForster Schreiber et al.|2011b[[Genzel et al.|2011[|Swinbank| 
|etal.20 1lfl 

A natural prediction of this model is that the old stars should 
be more smoothly distributed than the gas and young stars (some- 
what evident in Fig. [iJj^A related prediction is plainly that the 
clumps are relatively short-lived, with gas-phase lifetimes shown 
in Fig. [2} to some extent this should be manifest as young stellar 
population ages for the stars physically bound to the clumps. Test- 
ing these predictions requires observations that can probe the rest- 
frame near infrared structure of these galaxies. Moreover, patchy 
obscuration (given the very large gas densities) can make even 
the oldest background populations appear clumpy (see |Genel et al.| 
|2012| ). And as also noted in |Genel et aLl ( |2012| ), the fact that clumps 
form from both gas and stars means that - especially if the star for- 
mation efficiency is low and lifetimes are short - there will always 
be contributions to from pre-existing disk stars (with ^Gyr ages), 
which will introduce contaminating old stellar populations to all 
clumps as well as gradients in the clump ages with galactocentric 
radius (the clump age gradients simply tracing the pre-existing disk 
gradients). We therefore defer more detailed predictions of clump 
appearance in different observed wavelengths to future work, in 
which we will use radiative transfer models to construct mock 
observations. However, preliminary observational studies have al- 
ready shown that clumps dominate the light only in young stellar 
populations and represent transient local fluctuations in luminos- 
ity; when constructing stellar mass maps, clumps become far less 
significant ( |Wuyts et al.|2012 i A rigorous comparison of the sim- 
ulations here (using mock observations in the relevant observed 
bands) and these observations will be the subject of future work, 
but preliminary comparisons suggest that they strongly disfavor the 
"weak/no feedback" models (Wuyts et al., in prep). 

Which feedback mechanisms are most important in dissoci- 
ating clumps and recycling gas into the ISM? We have shown in 
Paper II and Paper III that in these dense, gas-rich systems radia- 
tion pressure is dominant over sources of thermal heating such as 
shocked SNe ejecta, stellar winds, and HII photoionization. This is 
because the gas densities within clumps are so large that the cooling 
time is significantly shorter than the dynamical time. This explains 
why previous simulations with feedback included only in the form 
of thermal heating from SNe have seen very little effect on the evo- 
lution of massive star-forming disks (see §[TJ. 

Within the context of radiation pressure, we also find that it is 
specifically the ability to drive winds out of clumps (and, possibly 
but not necessarily, the galaxy) that has the largest effect on long- 
term clump survival and bulge formation. This supports the conclu- 
sions of ' Genel et al.| ([20T2|, who found the same in fully cosmo- 
logical simulations, but with a phenomenological wind model (see 
§ [1]). Since the models here self-consistently produce winds with 

^ The more prominent rotation - since there is no vertical feedback support 
- in the no-feedback case leads to the mini-spiral structure seen in some 
cases in Fig.[T] for a more detailed discussion of the rotation support in the 
SNe-only case, we refer to |Ceverino et al.|j20T2| . 

This is always true to some extent, however it will plainly be moreso if 
the clumps themselves are short-lived and transient, rather than long-lived 
features. 



comparable global mass-loading factors (see Paper III), many of 
the conclusions should be the same (and we should expect to find 
similar results in cosmological simulations). In detail, the global 
wind mass-loading factors here are relatively modest, Mwind ~ 
(Fig. 6-7 in Paper III), a factor of ^ A lower than those in |Genel| 
|et al.| ( |2012| ) and closer to the values (at this mass and redshift) 



typically invoked in cosmological models to rep roduce the galaxy 
stellar mass function |Oppenheimer et al" 2010|p] But the implica- 



tions for clump survival are similar. The reason is that the feedback 
models here give rise to large local outflows and fountains (which 
do not appear in the |Genel et al.|2012l model), driving turbulence 
and outflows from individual clumps at modest velocities (several 
tens of kms~^). This is sufficient to unbind gas from a clump but 
not from the galaxy as a whole. Only a fraction of the material in 
these fountains is then further accelerated into the galactic outflow 
at > 100km s~^ A generic expectation of our more realistic mod- 
els is that there should be significantly more mass at modest ve- 
locities compared to the true "super-wind;" therefore even modest 
mass-loading factors of > 0.5M* at high velocities are consistent 
with our predictions here and probably correspond to sufficient lo- 
cal outflows and fountains to suppress clump inspiral. 

|Krumholz & DekeT|p010| ) used analytic models to consider 
whether or not radiation pressure would disrupt the clumps ob- 
served in high redshift galaxies, and concluded that it would not. 
There are a number of important differences between the full non- 
linear results here and their simplified model. The clouds here are 
clumpy and marginally bound, with a large fraction of the clump 
mass only transiently associated with a given cloud (as opposed 
to tightly bound or virialized; see Paper II). A combination of disk 
shear, gravitational instability and feedback operate simultaneously 
(lowering the threshold for efficient feedback). |Krumholz & Dekel| 
( |2010| ) assumed a minimal momentum flux from young stars (about 
the flux from the single scattering limit); considering the addi- 
tional contributions from SNe and stellar wind momentum flux, 
and multiple scatterings in the IR, the momentum flux here is a fac- 
tor of several larger. And they assumed a star formation efficiency 
of 1% per clump dynamical time, even though they argue that 
this fails to prevent collapse. In Paper I and Paper II, we show that 
whenever feedback is insufficient to disperse a GMC, runaway lo- 
cal collapse proceeds until the SFR increases to the level needed to 
once again further suppress collapse (see also ,Tasker|2011|[Dobbs| 
|et al.|20lT| ). Thus a higher efficiency always arises when feedback 
is "initially" insufficient. Overall, our calculations demonstrate that 
an instantaneous star formation efficiency of just ^ 1 — 5% (or an 
integral efficiency of ^ 5 — 20%) is more than sufficient to unbind 
(most of) the clumps. 

Given the stellar feedback mechanisms above, the models here 
predict local outflows from clumps in qualitative agreement with 
observations ( [Genzel et al.||201l"| [Newman et al.|[2b l2 ). Both our 
models here and those in [Genel etaL] ( |2012| ) indicate that, if the 
wind mass is comparable to that observed, clumps will, in fact, 
tend to be disrupted. This is also consistent with what is observed 
at low redshifts (albeit in different galaxy types): in some starburst 
and mergers the densities and velocity dispersions are similar to 
these high-redshift systems, and so the characteristic Toomre mass 
is similar; indeed, 100 GMCs with masses 10^ - 3 x lO'^ Mq 
and radii 0. 1 — 1 kpc have been observed in the local Universe in 



For a detailed comparison of the wind mass-loading, structure, veloci- 
ties, and phases with observations and cosmological constraints, we refer to 
Paper III. 
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such extreme systems (see 


Rand & Kulkarni|1990||Planesas et al.| 


|1991||Wilson et al.||2003| 


2006| and references therein). A com- 



bination of constraints from stellar population ages and masses, 
the maximum mass of star clusters observed in the same galaxies, 
the specific frequencies of bound/globular clusters, and the cluster 
formation rates all suggest that these clumps have (gas) lifetimes 
< lOOMyr, turn only a modest fraction ^ 10% of their mass into 
stars, and typically do not remain bound. In future work, we will 
extend our low-redshift models to starbursts and mergers to com- 
pare directly with these observations, but preliminary comparison 
suggests that they self-regulate in similar fashion. 

Our calculations do not include the inflow of gas via continu- 
ous accretion and mergers in a cosmological context. This allows us 
to robustly distinguish the effects of true isolated evolution and disk 
instability from e.g. ongoing minor mergers. But as a result, our 
conclusions should be viewed as statements about the relative effi- 
ciency of internal evolution with and without feedback, not state- 
ments about absolute bulge masses/structure. If the models here are 
representative of evolution over a few hundred Myr, then a minor 
merger of mass ratio 10: 1 or larger would be able to induce compa- 
rable bulge growth, and a major merger would change the system's 
properties by more than several Gyr of evolution in this mode. Con- 
tinuous accretion may, however, accelerate the evolution produced 
by the short-lived clumps and/or effect angular momentum transfer 
(see |Keres et al.|200 5^). In future work, we will extend these models 
to cosmological zoom-in simulations to explore these questions in 
more detail. 
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APPENDIX A: NUMERICAL EFFECTS OF THE 
RADIATION PRESSURE IMPLEMENTATION 

The numerical details and tests of the feedback models are dis- 
cussed in great detail in Paper I, Paper II, and Paper III. However, 
since we argue that radiation pressure is the most important element 
of feedback in these dense systems, we here discuss how our main 
conclusions change if we vary numerical aspects of the radiation 
pressure model implementation. 

In Pig. |A1| we plot the final mass profiles in the simulations 
as in Pig. |3] but vary the radiation pressure physics included. The 
"standard" (all feedback mechanisms included) case and "no radi- 
ation pressure" case are the same as Pig. [3] 

Recall, the radiation pressure force is coupled in two sepa- 
rate algorithms, the "local IR" and "long-range" terms described 
in §[2] (mechanisms [1] and [5], respectively). The former includes 
the effect of primary optical/UV photons absorbed in the immedi- 
ate vicinity of young stars (within of order a smoothing length), as 
well as multiple- scattering of the re-emitted IR photons when the 
IR optical depth r is significant (which dominate this term when 
r > 1). The latter includes the effect of the photons escaping this 
region (UV/optical, and IR) as they propagate to infinity; although 
the IR term is included, since most of the volume outside the dense 
regions near stars is optically thin to it, most of the "work" in this 
term is done by the UV/optical photons being absorbed by less- 
dense gas further from the dense regions of origin. It is instructive 
to decompose the effects of these terms. 

In the Appendix of Paper II, we discuss at length the accu- 
racy of our long-range radiation pressure approximation. There, we 
show that although we are not solving the full radiation hydrody- 
namic equations, the simplifications used here are a surprisingly 
good approximation to full radiative transfer solutions. Neverthe- 
less, it is important to check how big an effect this has. In Fig. |Al| 
we therefore re -run our simulation, but instead of using the on-the- 
fly simulation information on stellar populations and gas distribu- 
tions to calculate the spectral shape escaping and propagating in the 
long-range radiation pressure term, we simply fix the SED shape to 
one chosen to match observations. The details of this procedure 
are given in the Appendix of Paper II; as noted there, choosing in- 
stead a mean SED matched to a full radiative transfer calculation 
on the simulations gives an almost identical result (since it closely 
matches what is observed). This has almost no effect on our con- 
clusions relative to our "standard" model. 

If we remove the long-range terms (mechanism [5] in §[2| en- 
tirely, but still include the local/IR term (mechanism [1]), we see 
a much more significant difference. Both the bulge mass and in- 
ternal evolution of the outer disk are significantly enhanced above 
the "standard" model. Most of the new bulge mass is at ^ 1 kpc. 
The very central regions at < lOOpc are regulated by the IR term 
and so are still maintained at low densities, but the overall ability to 
disrupt the kpc- scale clump complexes is reduced. The net effect 
is comparable to - although still clearly less dramatic than - the 
effects of removing radiation pressure entirely. 

It therefore seems that the long-range radiation acceleration is 
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at least comparable, if not more important, to the local IR radiation 
pressure in suppressing the clump star formation efficiencies and 
subsequent bulge growth. This is true for several reasons: but the 
critical point is that the kpc-scale complexes have clear substruc- 
ture. As shown in Paper II, at any given instant, most of the mass in 
complexes is not in the extremely dense sub-regions forming stars, 
and this "intermediate-density" material is essentially never tightly 
bound or virialized inside a cloud when feedback is present. The IR 
radiation pressure, obviously, acts significantly only where the IR 
opacities are substantial. The volume- average surface densities of 
the complexes are only marginally optically thick to the IR (tir ^ 1 
for Egas lO^M0kpc~^); so the IR coupling is negligible in the 
lower-density volume-filling medium inside each complex. The IR 
coupling is therefore dominated (as shown in Paper I) by very 
dense sub-clumps at 1 pc scales (densities Sgas ^ 10"^ Mq kpc~^ 
and three-dimensional densities n > lO^cm"^), over a very short 
time when they have just formed young stars. But these dense sub- 
regions are distributed throughout the complex, so the IR term does 
not act coherently across the entire cloud complex to disrupt it. 
Rather, it disrupts the dense star-forming sub-regions shortly after 
they form stars, recycling the dense gas into the parent complex, 
at turbulent velocities of order the local escape velocity from the 
sub-clump ('-^ 10 — 50kms~^). This makes it optically thin to the 
IR photons, but (being marginally bound) much more vulnerable to 
continuous, coherent acceleration from the collective UV/optical 
flux of the all the stars in the complex. 

The short-range radiation pressure mechanism - in particular 
the multiple- scattering of IR photons producing a local momentum 
flux of ^tL/c- can, in principle, be more affected by the details of 
accurate radiation transport. For example, [Krumholz & Thompson] 
( |2012| ) argue that, in a time-averaged sense, the "boost factor" in 
p^TeL/c does not scale linearly with the full average optical depth 
T when T ^ 1 and does not typically exceed a factor of several. We 
should stress, however, that our conclusions are actually not that 
different from those of [Krumholz & Thompson] ( |2012| ). Those au- 
thors argue specifically that the reason is "capped" is because, 
whenever the radiation pressure force temporarily exceeds gravity, 
gas is locally blown away and dense regions are "shredded" by ra- 
diation Rayleigh-Taylor instabilities, until the coupling is lowered 
such that radiation balances gravity. Even though we do not fol- 
low the full radiation transport, this is the generic behavior in our 
simulations as well (IR radiation pressure drives to lower effective 
Te to offset gravity). As a result, as discussed in Paper I and Pa- 
per II, the IR term primarily drives local turbulence, disrupting or 
shredding dense sub-regions and driving turbulence to maintain a 
Toomre 2 ~ 1, while the long-range radiation pressure (which op- 
erates via single scattering and is not affected by this argument) 
drives winds in loosely-bound material, primarily coupling via its 
UV/optical opacity. Consider, for example, in Paper I, we show that 
the mean true IR r in the dense protostellar core-like regions that 
have just formed new stars - which dominate the momentum input 
- can reach as high as ^ 10 — 50. However, we also show that the 
coupling of tL/c in these regions immediately (often in < IMyr) 
disperses the gas to lower r around the source and allows resolved 
leakage, so that the luminosity -averaged coupling factor, i.e. the 
true Te we should compare to [Krumholz & Thompson] ( |2012| ), is 
actually just 4 — 5 in our simulations. Both this and the result- 
ing local velocity dispersions and Eddington ratios are in reason- 
able agreement with these and other radiative transfer calculations 
( [Kuiper et al.|2012| Jons son et al. in preparation). 

Still, one could imagine sub-resolution scale effects that might 
reduce this local IR coupling. We discuss this extensively in Pa- 



per I (Appendix B), and in particular consider the case of a fractal 
ISM continuing to infinitely small scales below the resolved region 
used to calculate r; as opposed to a homogenous medium, this will 
present "holes" or "channels" along which photons can leak out. 
For a lognormal or power-law density distribution, we show that 
in this particular case it is straightforward to analytically calcu- 
late the "effective" coupling and leakage effects. In all cases where 
r is large, this dramatically enhances the escape fraction of UV 
and optical photons. However, if the scatter in density contributed 
strictly below the resolution scale of ^ 1 — lOpc is weak or mod- 
est, < 0.5 dex, then this has no effect on the dimensional scaling of 
the momentum imparted by multiply- scattered IR photons (rather it 
just amounts to a weak normalization correction). For larger scatter 
it introduces a scaling Te ^ t^ where /3 < 1 is a weakly decreas- 
ing power of the sub-resolution density dispersion. In Paper I we 
show how this can be implemented in the simulations and used to 
allow both for the corrected IR coupling and UV/optical escape 
fractions. We therefore consider two examples using this "leakage" 
estimator. In the first, we follow Paper I (model "PL: a calculated" 
therein) and use the local (resolved) density dispersion added in 
quadrature with a large assumed floor of ^ 1 dex and the solution 
for a fractal medium - this results in a highly sub-linear re(r) (ap- 
proximately oc T^/^), and actually produces significantly weaker Te 
values than those calculated in full radiative transfer calculations by 
either ^umholz & Thompson] pOTI] ! or lKuiper et"aLl ( [20T2l ). But 
we see little difference between this and our standard model - if 
anything, the bulge is slightly smaller! This occurs because, as dis- 
cussed above, the single- scattering is critical to actually dispersing 
the cloud complexes, while the IR largely drives local turbulence 
and is always driven to approximate equilibrium with gravity. The 
latter equilibration still occurs, even if the coupling is made weaker 
(gas collapses to slightly higher densities first, giving higher r to 
obtain the same r^); this is shown explicitly (along with the lack of 
strong effects on any of our predicted galaxy properties) in Paper I 
(Figs. 5,9, 10, B2). And allowing for leakage of course makes the 
long-range radiation pressure forces stronger, by removing fewer 
of the UV/optical photons in the dense regions near stars where 
they act less efficiently. 

To clearly distinguish between the effects of decreasing the 
local IR radiation force and increasing the long-range force, we can 
artificially multiply the local IR radiation pressure by an arbitrary 
factor, without increasing the UV/optical leakage. Fig. ]A1] shows 
such an experiment, with a reduction of the local IR force by a 
factor of 3, which leads to a very modest galaxy-averaged Te ~ 2. 
As Fig. ]A11 shows, here we do find more bulge growth and longer- 
lived clouds than in the previous case, as expected, although the 
effect is small. 

Other effects discussed in Paper I, however, could mean that 
our implementation actually under-estimates the strength of the ra- 
diation pressure terms. For example, if the gas density profiles rise 
sufficiently steeply around young stars, or if we are under-resolving 
the maximum densities (hence opacities) to which regions sur- 
rounding very young stars collapse, then our simulation-calculated 
momentum flux could be a substantial und er-estimate. | Kuiper et al. 
(^012 ) also argue that calculations such as ]Krumholz & Thompson 
( [2012^ , which do not follow the exact radiative transfer but adopt a 
pure flux limited diffusion approximation and see efficient dispersal 
as described above, can underestimate the true coupling efficiency 
by an order-of-magnitude. Unsurprisingly, we see further reduction 
of bulge growth and internal evolution in experiments where we en- 
hance the radiation pressure strength to match these suggestions. 

Finally, as shown in Paper I and Paper II, we have run reso- 



© 0000 RAS, MNRAS 000, 000-000 



Feedback & Clump Coalescence 1 1 



12 




0.1 1.0 10 100 

R [kpc] 

Figure Al. Remnant stellar mass profiles as Fig. |3] here we specifically 
compare models where we freely vary sub-elements of the radiation pres- 
sure implementation. The "standard/all feedback" model and "no radiation 
pressure model" are identical to Fig. [3] The "long-range SED fixed to obs." 
model fixes the spectral shape of the propagated long-range radiation pres- 
sure (RP) terms (mechanism [5] in §[2| to the mean observed in the galaxies 
modeled. This has little effect, since the self-consistently calculated SEDs 
agree well with the observations The "no long-range RP" model removes 
the long-range terms entirely. This significantly suppresses the mixing and 
disruption of the largest super-clump complexes, resulting in more bulge 
growth and internal evolution. The "local IR RP force xO.3" model sup- 
presses the local IR RP term (mechanism [1] in §|2j by an arbitrary factor 
= 0.3. This leads to less feedback-regulation in the dense gas and greater 
bulge growth, but the effect is weak when the long-range RP is still present. 
The "local IR RP with strong leakage" model calculates the local IR RP 
term assuming a sub-grid fractal ISM with a broad density PDF giving 
a much larger UV/optical escape fraction and weaker multiple-scattering 
"boost" for a given average dust optical depth in the region. Although this 
makes the local RP term weaker, the enhanced escape of UV/optical pho- 
tons from the very dense regions around young stars makes the long-range 
RP forces significantly stronger and increases the coherent acceleration of 
gas between the dense sub-cloud cores in the massive cloud complexes (as 
well as the launching of RP-driven winds); as a result this slightly sup- 
presses the internal evolution relative to the standard model. 

lution tests with these models up to ^ 10^ particles and find very 
good convergence at the standard resolution adopted here; in fact 
for the bulk properties here we see reasonable convergence even at 
order-of -magnitude lower particle numbers. 
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